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ABSTRACT 

We show how the observed gamma ray burst (GRB) optical afterglow (OA) and red- 
shift distributions are changing in time from selection effects. For a subset of Swift 
triggered long duration bursts, we show that the mean time taken to acquire spec- 
troscopic redshifts for a GRB OA has evolved to shorter times. We identify a strong 
correlation between the mean time taken to acquire a spectroscopic redshift and the 
measured redshift. This correlation reveals that shorter response times favour smaller 
redshift bursts. This is compelling evidence for a selection effect that biases longer 
response times with relatively brighter high redshift bursts. Conversely, for shorter 
response times, optically fainter bursts that are relatively closer are bright enough 
for spectroscopic redshifts to be acquired. This selection effect could explain why the 
average redshift, < z >« 2.8 measured in 2005, has evolved to < z >xs 2, by mid 
2008. Understanding these selection effects provides an important tool for separat- 
ing the contributions of intrinsically faint bursts, those obscured by host galaxy dust 
and bursts not seen in the optical because their OAs are observed at late times. The 
study highlights the importance of rapid response telescopes capable of spectroscopy, 
and identifies a new redshift selection effect that has not been considered previously, 
namely the response time to measure the redshift. 
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1 INTRODUCTION 

GRBsQ are extremely bright transients observed out to great 
cosmological distances. Because of their high luminosity in 
7-rays, they are a unique probe to their host galaxies in the 
high-z Universe. At least a fraction of GRBs have been asso- 
ciated with the collapse of massive stars via the association 
of supernova signat ures observed with the fading GRB op - 
tical afterglow e.g. (|Hiorth et all 120031 : Istanek et al J [2003V 
The afterglow most likely originates from an external shock 
produced as the blast wave from the GRB collides with the 
interstellar medium (ISM). The successive afterglows at pro- 
gressively longer wavelengths (X-ray, optical, radio) are a re- 
sult of the expanding shock wave interacting with the ISM, 
causing it to slow down and lose energy. The GRB-supernova 
connection implies that GRBs should track the star forma- 
tion rate (SFR) of massive stars, and could be used as a 
complementary probe of the SFR in the high-z regime where 
optical data is lacking. 

With the launch of NASAs Swift satellite in 2004 



November, a new era of rapid GRB localization was born. 
A slew of detectors on board Swift provided the means to 
rapidly find GRBs with small error boxes, enabling rapid 
follow up of the optical afterglow by dedicated ground- 
based telescopes. Prior to Swift, only about 50% of local- 
ized GRBs were identified with an optical afterglow. The 
high sensitivity of Swift coupled with the growing number 
of rapid response ground-based telescopes capable of spec- 
troscopy promised to fill the gaps. Surprisingly, this did not 
happen: optical/NIR afterglows have been found for nearly 
80% of GRBs but only 4 0-50% have measured redshifts 
l|Tanvir fc Jakobssorj|2007h . It is now established that some 
GRBs may not have an afterglow and are defined as 'dark' 
bursts. Jakobsson et al. (2004) define a dark burst by sug- 
gesting that the spectral slope between the optical and X-ray 
is /3ox < 0.5. Fynbo et al. (2008) point out that up to 25% of 
Swift bursts fulfill this condition-see also lRol et all (2005). 



2 DARK GRBS AND SELECTION EFFECTS 
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We are now confronted with another problem - the i s sue o f 
the missing GRB optical afterglows. iLe Floc'h et all (|2006l ) 
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show out that Spitzer observations of long GRB are not asso- 
ciated with the massive and luminous infrared galaxies that 
are associated with the bulk of star formation in the early 
Universe. This work supports claims arguing for a GRB host 
population dominated by blue, young and low-mass galaxies 
and that dust extinction would not play a prominent role in 
such GRB hosts. However, there is some evidence that dust 
extinction in the local ISM of the host along the line of 
sight is the cause of some of the dark bursts. Evidence for 
this is shown in the study of the host galaxy of GRB 030115 
|Levan et al.ll2006h . Its optical afterglow was fainter than 
many upper limits for other bursts, suggesting that without 
early NIR observations it would have been classified as a 
dark burst. Both the colour and optical magnitude of the 
afterglow probably arise from dust extinction and indicate 
that at least some optical afterglows are very faint due to 
dust along the line of sight. We must conclude that at least 
some fraction of dark bursts are caused by dust extinction 
in the host galaxy . 

Roming et al. (2006), used very early observations of 
Swift GRBs to investigate the most probable cause of opti- 
cally dark GRBs. They find that ~ 25% of the afterglows 
in their sample are extincted by Galactic dust, ~ 25% are 
obscured by absorption in the immediate GRB environment 
and ~ 30% are most likely attributable to Ly-alpha blan- 
keting and absorption at high redshift. So it is highly likely 
that the dark bursts result from a combination of extinction 
factors relating to t he GRB environment, host galaxy type 
and galaxy redshift dMirabel fc Rodriquesll20"03h . Schady et. 
al. (2007) show that many OAs have significant optical ex- 
tinctions compared to the optically bright bursts and argue 
that extinction could account for many of the OAs not ob- 
served by the Ultr a- Violet Optical Te lescope onboard Swift. 
In support of this, iTanvir et alj (|2008l ) show that dust is the 
likely cause of the colou r and faintness of the red OA of GRB 
060923A. Furthermore. I Jaunsen et al] {2008) show that the 
highly extinguished afterglow of GRB 070306 is due to host 
galaxy dust extinction. 

The probability of obtaining a reliable GRB redshift is 
directly related to the signal-to-noise ratio of the absorption 
or emission lines. Ideally, multiple 'strong' lines are required, 
but this is made difficult because GRB OA brightness de- 
cays rapidly in time. Combined with the fact that many 
GRB host galaxies are too faint for redshifts to be obtained, 
the time taken to image the OA with medium to large tele- 
sc opes become s a cri tical factor. This was first pointed out 
bv lFiore et ail \200fl ) for the observed discrepancy between 
the HETE and BeppoSAX redshift distributions compared 
to Swift. 

It is clear that these optical afterglow selection effects 
will propagate through to the observed redshift distribu- 
tion. In addition to afterglow extinction biases above, the 
so-called 'redshift desert' in z ~ 1 — 2 is a region where 
it is difficult to meas ure redshifts bec ause of the lack of 
strong emission lines (|Fiore et al]|2007f ). However, we note 
that for absorption lines, the Mg II doublet is prominent in 
0.4 < z < 2.2, so that the redshift desert may not play such 

a prominent role in redshift determination. 

These problems were furth er highlighted by ICow ard 

l|2007h and ICoward et al] (|2008h . They argued that in z = 
— 1, the GRB redshift distribution should increase rapidly 
because of increasing differential volume sizes and strong 



SFR evolution. Until mid 2007, this characteristic in the 
Swift redshift distribution was not apparent. To account for 
this discrepancy, they argue that other biases, independent 
of the Swift sensitivity, are required. The lac k of measured 
redshi fts in z £s 1 — 2 up to 2007, discussed bv lCoward et al] 
(2008), is surprising and has not been resolved. 

2.1 GRB redshifts and the 'Learning Curve' 

Assuming a redshift can be measured from any burst in the 
Universe and that GRBs are independent events (see How- 
ell et al. 2007), then the bursts will follow a Poisson distri- 
bution defined by their mean rate observed in our frame. 
Furthermore, the mean GRB redshift will fluctuate about a 
constant assuming an observation duration much less than 
any cosmological epoch. This also implies that the observed 
time-series comprised of GRB redshifts will be stationary in 
time i.e. the statistical moments will not drift in time. Small 
number statistics will obviously result in scatter of the sta- 
tistical moments, but cannot be the cause of any evolution 
in the statistical moments of the GRB redshift distribution 
on a time-scale of years. 

A shift of the mean of the GRB redshift dis tribution was 
obser ved in the early part of the Swift mission (|Berger et al] 
2005). This was attributed to the improved sensitivity and 
more accurate localisation by Swif t resulting in a b i as for 
fainter and higher redshift bursts. Ijakobsson et al] (|2006T ) 
showed that within the first year of Swift the mean redshift 
for a subset of 28 bursts had drifted to about 2.8, about 
double that of the pre- Swift average redshift. 

Assuming that satellite sensitivity is the dominant fac- 
tor for determining redshift statistics, one could assume that 
the high mean redshift observed in the early part of the mis- 
sion would remain fairly constant. If other factors impact on 
the statistics, such as the time taken to acquire high signal- 
to-noise spectra, the statistics may well reflect this. Given 
that Swifts sensitivity and GRB localisation ability has not 
degraded over time, one must consider the next link in the 
chain for measuring GRB redshifts: optical follow up by tele- 
scopes capable of spectroscopy. 

Given the above argument, any gross time-dependent 
evolution observed in the mean of the GRB redshift distri- 
bution must result from how the redshifts are measured. As 
discussed above, this is not a straightforward process, but 
depends on obtaining high quality spectra from OAs that 
are decaying rapidly in time. We will show in this Letter 
that the time taken to measure a spectroscopic redshift can 
cause an important selection effect. By analysing the delay 
time for spectroscopic redshift data as time-series, we show 
how the non-stationarity of the data can be used to con- 
strain the origin of biases that plague this data. This type 
of analysis can help discriminate between ground-based op- 
tical selection effects and possible intrinsic characteristics of 
GRBs. 

The learning curve effect, first mathematically modelled 
bv lWrightl (|1936h . has a strong influence in any large organi- 
zation that is technology and skill dependent. It is based on 
the idea that efficiency for a particular task improves with 
experience and technology. It has become an important tool 
in management to model how productivity and effici ency 
improves in large organisations l|Adler and Clarklll99lf ). 

With the globalization of optical astronomical tech- 
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niques and knowledge, the learning curve effect is expected 
to influence the efficiency of targeted transient surveys, such 
as supernova search programs (Coward 2008 in preparation). 
For the case of GRB OA spectroscopic redshift measure- 
ments, the effect is very subtle but nonetheless important. 
This basic idea can be related to GRB optical follow-up 
in terms of the efficiency of localizing, imaging and obtain- 
ing high quality spectra suitable for measuring redshift. In 
particular, a critical factor that determines the efficiency 
of these tasks is the time it takes to optically localize the 
rapidly fading OA and to obtain high signal-to-noise spec- 
troscopy. One would expect that this task should be improv- 
ing in time and we seek to identify if this effect is impacting 
on the redshift distribution statistics as a function of time. 



3 DATA ANALYSIS AND RESULTS 

We select 192 GRBs detected by Swift from 2005 March to 
2008 July with high-energy emission duration greater than 4 
s. After removing those bursts without redshifts, the remain- 
ing 91 consist of those with absorption and emission spectra 
of the OA. From this population, we select a subset of 64 
from GCN circulars that have response times from when the 
burst was triggered by Swift's BAT to the acquisition of a 
spectroscopic redshift. We have attempted to excluded burst 
redshifts measured from the host galaxy at very late times 
as obtaining these spectra did not depend critically on the 
response time of the telescope. 

We employ two additional statistics for the analysis: the 
optical brightness at discovery and the time of the burst. 
The optical brightness is not measured with the same fil- 
ters across the data. Nonetheless it is a useful parameter for 
investigating if the efficiency of localizing bursts by rapid 
follow up of optical telescopes is changing in time. The red- 
shift data above are analysed as time series to probe how the 
statistical moments-i.e. the mean, variance and the discov- 
ery rate-evolve over the mission time of Swift. To determine 
if the data is non-stationary in time, a moving average filter 
spanning 4 nearest neighbours on each side of an event is 
employed. 

Fig. 1 plots the raw redshifts and output from the mov- 
ing average filter (using 4 nearest neighbours) as a time- 
series. The associated GRB OA magnitudes at discovery 
are also plotted. Although there are fluctuations over pe- 
riods of months, there is a clear downward trend in the red- 
shift averages over a 3 year period. We test if the observed 
non-stationarity of the redshift time-series is related to how 
the redshifts are measured, in particular the time taken to 
obtain spectroscopic redshifts. 

Fig. 2 plots both the response times to acquire these 
64 redshifts and a moving average of this data, < T z >, 
against the time when the burst occurred. The plots shows 
a definite long term trend in the response times. To test if 
< T z > is affecting the selection of GRBs in a certain dis- 
tance range, Fig. 3 plots < T z > against the mean redshift. 
The plot provides compelling evidence that there is a signif- 
icant selection effect at work. There appears to be a strong 
correlation between the response times to acquire a redshift 
and the redshift measured via spectroscopy of the OA. 

Interpreting this result provides the key to understand- 
ing the observed drift in the mean GRB redshift over time. 



It is clear that the probability of observing smaller GRB red- 
shifts has increased up to mid-2008. From Fig. 3, the shortest 
waiting times to acquire a redshift correspond, on average, to 
smaller redshift bursts. This observation can be understood 
if these bursts have rapidly fading OAs that would have been 
too faint for high quality spectroscopy if observed at later 
times. Conversely, the longer duration response times prefer- 
entially select those bursts with very bright OAs that can be 
observed spectroscopically at higher redshift. Furthermore, 
this selection effect can be reconciled with the previously ob- 
served lack of redshifts in z < 2, as noted by several studies 
including Coward et al. (2008). 

Fig. 4 plots the cumulative growth in measured redshifts 
as a function of observation time for the redshift regimes z > 
2 and z < 2. The plots show that the number of high- 2 bursts 
measured is roughly constant, compared to the marginally 
non-linear increasing rate of the z < 2 events. It is apparent 
that the probability of measuring a z < 2 burst is increasing 
with time, a result supporting the observed drift to a smaller 
mean redshift shown in Figures 1 and 3. The root of this 
selection effect may originate from intrinsic differences in 
the OA brightness (at the source) between the low-medium 
and high redshift bursts. 



4 SUMMARY AND DISCUSSION 

The results based on the analysis are summarized below: 

1. A time-dependent trend in the mean of the redshift dis- 
tribution is observed. The mean of redshift has drifted from 
about 2.8 in 2005 to about 2.1 by mid-2008 (see figure 1). 

2. Over the same observation period we show that the mean 
time taken to acquire spectroscopic redshifts from the GRB 
OA has become shorter. In the period 2005 to 2006 Nov we 
find that the average time to localize and acquire a spectro- 
scopic redshift was about 800 min. In contrast, the response 
time in 2007 July to 2008 July has reduced to about 200 
min. In 2008 large telescopes such as the VLT have man- 
aged to acquire OA localization and spectroscopic redshifts 
in under 100 min for half of the bursts observed. 

3. There is a clear correlation between the mean time taken 
to acquire a spectroscopic redshift and the measured red- 
shift. Specifically, shorter response times correspond to, on 
average, smaller redshifts (see Fig. 3). This is compelling ev- 
idence for a selection effect that biases optically bright high 
redshift bursts with longer response times. Conversely, for 
shorter response times, optically fainter bursts that are rela- 
tively closer will be bright enough for spectroscopic redshifts 
to be acquired. 

4. The improved efficiency in the response times to acquire 
redshifts is evidence that the learning curve effect is having 
an impact on the redshift statistics. Interestingly, the effect 
is not manifesting as a significant growth in the total num- 
ber of measured redshifts over time. It does appear to be 
introducing a time-evolving bias for selecting OA at smaller 
redshifts (see Fig. 4) compared to earlier Swift detections 
prior to 2006. 

It is interesting that the apparent deficit of z < 1.5 
bursts, evident in the redshift distribution prior to 2007 e.g. 
Coward et al. (2008), is less pronounced in 2007-2008. This 
could be a direct result of bursts in z < 1.5 having relatively 
fainter OAs compared to higher redshift bursts. As the tele- 
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Figure 1. The subset of 64 GRB redshifts with confirmed re- 
sponse times for spectroscopy (red squares) plotted as a time- 
series, starting in March 2005 and ending mid-2008. Output from 
a moving average filter for the same data is plotted (solid line), 
showing how the average redshift has a time-dependant trend to- 
wards decreasing values. The observed trend cannot be a result of 
evolution of the sources, but must be related to how the redshifts 
are determined. Optical magnitudes at discovery (blue circles) 
are also plotted, revealing a weak trend for observing increasingly 
brighter bursts at early times. 



10' 



10 



10' 




06 



07 
year 



08 



Figure 2. Plot showing the evolving reduction in average re- 
sponse times, < T z >, to acquire a spectroscopic redshift from 
the GRB OA as a function of Swift mission time (solid line). The 
raw data (solid squares), taken from GCN circulars, consists of 
64 response times measured from the Swift's BAT trigger to the 
time of obtaining a high quality spectrum. The data span from 
2005 March to mid 2008. Average response times have evolved 
from about 1000 min in 2005 down to several hundred minutes 
by mid-2008. This increase in efficiency is an indicator of the 
'learning curve' effect at work. 



scope response times have reduced, there are increasingly 
more fainter OA that are bright enough for spectroscopic 
measurement. We are currently investigating if the drift in 
the redshift distribution statistics can be used to probe the 
underlying reason why z < 1.5 bursts have been relatively 
more difficult to determine. In our analysis, we have not at- 
tempted to construct a model to account for the shifting of 
average redshift measured with response time. In a follow- 
on study, we plan to investigate the correlations found in 
this study via simulation and by fitting a model to the aver- 
aged redshift-response time data. This will provide a means 
to confirm the claim in this study that GRBs at high- z are 
intrinsically optically brighter relative to those bursts at low- 
medium redshift. 

Finally, we highlight the critical role that telescopes ca- 
pable of rapid response spectroscopy have on probing the en- 
vironments of GRBs over different cosmological epochs. This 
work shows that to obtain an unbiased sample of redshifts, 
will require understanding a new bias that has not been con- 
sidered previously, namely the response time to measure the 
redshift. 
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Figure 3. Plot showing the correlation of average response times, 
< T z >, to acquire a spectroscopic redshift with the mean of the 
redshift distribution. The plot shows that the shortest response 
times on average correspond to smaller redshift bursts. This im- 
plies that these bursts have rapidly fading OAs that would have 
been too faint to acquire high quality spectroscopy from if ob- 
served at later times. Conversely, the longer duration response 
times may cause a preferential selection of high- 2 bursts that have 
very bright OAs and miss smaller- z bursts if they have relatively 
fainter OAs. For < T z >Ri 1000 min, the average redshift is about 
3, and for < T z >ss 200 min, corresponding to observations in 
2008 (see Fig. 2), the average redshift has reduced to about 2. 
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Figure 4. Plot showing the cumulative growth in measured 
redshifts over the time of the Swift mission for the two redshift 
regimes: z < 2 (solid line) and z > 2 (dashed line). The growth in 
z > 2 redshifts is approximately linear while the numbers of z < 2 
redshifts is increasing at a marginally non-linear rate, mainly be- 
cause of an enhancement in 2008. 
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